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The vortex matter in bulk type-II superconductors serves as a prototype system 
for studying the random pinning problem in condensed matter physics. Since the 
vortex lattice is embedded in an atomic lattice, small angle neutron scattering 
(SANS) is the only technique that allows for direct structural studies. In tradi- 
tional SANS methods, the scattering intensity is a measure of the structure factor 
averaged over the entire sample. Recent studies in vortex physics have shown that 
it is highly desirable to develop a SANS technique which is capable of resolving 
the spatial inhomogeneities in the bulk vortex state. Here we report a novel slic- 
ing neutron diffraction technique using atypical collimation and an areal detector 
which allows for observing the three dimensional (3D) disorder of the vortex mat- 
ter inside an as-grown Nb single crystal. 



A major goal in the field of vortex physics is to determine 
the nature of the long-range order in the ground state of the 
Abrikosov vortex line lattice (Abrikosov, 1957). It has been the- 
pretically predicted that in systems with weak random pinning, 
a new phase of matter known as a Bragg glass exists in the 
vortex state (Chudnovsky, 1989; Nattermann, 1990; Giamarchi 
& Le Doussal, 1995). In contrast to a perfect crystal with a 
delta-function structure factor, the Bragg glass theory predicts 
a power-law shaped structure factor (Chudnovsky, 1989; Gia- 
[marchi & Le Doussal, 1995; Bogner et ai, 2001), similar to a 
'2D crystal at finite temperatures with no random pinning (Imry 
,& Gunther, 1971; Nelson & Halperin, 1979). Since the vortex 
•lattice is embedded in an atomic lattice, neutron scattering (Fur- 
jer et ai, 2009) is the only technique for direct structural stud- 
ies. There have been several attempts (Klein et ai, 2001; Dani- 
ilidis et ai, 2007; Laver et ai, 2008) to observe the Bragg glass 
phase, but no conclusive evidence has been found even though 
there are compelUng arguments (Klein et ai, 2001) for why one 
would not be able to measure it directly. 

In particular, the vortex states in bulk Nb samples have been 
studied extensively (Schelten et ai, 1971; Christen et ai, 1975; 
Lynn et ai, 1994; Rosov et al, 1994; Ling et ai, 2001; Forgan 
et ai, 2002; DaniiHdis et ai, 2007; Laver et al., 2008). These 
SANS measurements revealed the bulk properties of the interior 
of the samples. However, it has been recently observed that the 
vortex matter in a Nb single crystal can be spatially inhomo- 
geneous (Hanson et ai, 2010). Since the SANS measurement 
is based on an averaged signal, it would be difficult to discern 
the nature of disorder in the vortex matter It is thus necessary 
to develop a tool to spatially resolve the vortex matter structure 
from different locations in a bulk superconductor. 

By using a thin ribbon-shaped beam, as prepared on the 
Advanced Neutron Diffractometer/Reflectometer (AND/R) at 



the NIST Center for Neutron Research (NCNR), we find that 
it is possible to spatially resolve distinct, locally coherent 
regions of the vortex lattice within the sample (similar to crys- 
tallographic grains, with a small angular mismatch between 
regions), transverse to the direction of the magnetic field (Han- 
son et ai, 2010). Here, we report an improvement on this 
method which allows us to probe variations of the the in-plane 
vortex matter structure factor along the direction of the applied 
magnetic field. We demonstrate the effectiveness of this method 
in reveaUng the splitting of a Bragg peak along the magnetic 
field direction of the vortex lattice in a Nb single crystal. 

2. Sample Properties 

The sample is an as-grown Nb single crystal (cylinder shaped) 
with its (111) direction along the cylindrical axis, used in the 
previous study (Hanson et ai, 2010). It has a diameter of 
12.1 mm, a height of 10.1 mm, and the weight is 9.69 g. The 
residual resistivity ratio (RRR) is 32. By using the ac magnetic 
susceptometry technique (Ling & Budnick, 1991), we have 
mapped out the phase diagram and observed a weak peak effect 
in this sample. The characteristic field values at T = 4.2 K are 
the lower critical field Hd — 0.100 T, the upper critical field 
Hc2 = 0.490 T, and the field at the peak of the peak effect 
Hp = 0.400 T. 

We apply the magnetic field parallel to the (111) direction of 
the Nb crystal. The vortex lattice is prepared with a zero-field 
cooling procedure followed by thermal annealing and a small ac 
magnetic field perturbation. The magnetic field {H = 0.140 T) 
is applied after the sample is cooled through the transition 
temperature of 7^. = 9.2 K down to T — 4.0 K. We thermally 
anneal the vortex matter by slowly increasing the tempera- 
ture (in increments of AT — 0.5 K) to T = 7.5 K, just below 
the transition temperature, Tc2(H — 0.140 T), and then slowly 
cooling back down to T = 4.0 K (here AT" = 1.0 K). Finally, 
we perturb the vortex matter by applying an ac magnetic field 
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(sawtooth-function) with an amplitude of 0.001 T and a fre- 
quency of 1.0 Hz for 5.0 minutes provided by a Helmholtz coil 
outside the superconducting magnet. We find that this proce- 
dure is optimal for forming a stable vortex state with the highest 
order in this sample. 

3. Conventional Neutron Scattering on Vortex Matter 

In a neutron scattering experiment, the Bragg condition (Bragg 
& Bragg, 1913) occurs when the scattering vector Q equals the 
reciprocal lattice vector G, where \G\ = ^ = ^ sin 9. Here A 
is the neutron wavelength, 6 is the Bragg angle, and d is the 
distance between lattice planes. In vortex physics (Tinkham, 
1996), the apphed magnetic field, H, determines the value of 
d through the relationship 

where a is the vortex lattice constant and 3>o = 2.07 x 10~^^ Wb 
is the flux quantum. The sketch of an ideal vortex lattice in our 
sample is shown Fig. 1(a). 

The AND/R scattering geometry is shown in Fig. 1. The 
parameters of the setup are described elsewhere (Dura et ai, 
2006) but some important aspects are reiterated here. There are 
two slits in front of the sample, each one defined by two borated 
aluminum plates. With both slit widths set to 1.0 mm, the inci- 
dent neutron beam has an angular spread of HWHM= 0.042°. 
A wide third sUt is placed directly behind the sample to reduce 
the background scattering that reaches the areal detector. A 
fourth slit can be placed directly in front of the detector to fur- 
ther reduce background noise and improve the instrument reso- 
lution. In Fig. 1 (b), only the second slit is plotted. 

The azimuthal direction (perpendicular to G) is probed by 
rotating the sample and the magnetic field (and the associated 
reciprocal lattice) while keeping the detector (and the associ- 
ated scattering vector) fixed in space at an angle 20. The radial 
direction (along G) is measured by rotating the sample (and the 
magnetic field) and neutron detector in concert, always moving 
the detector by 2Aw for every sample movement of Aw. This 
procedure keeps the scattering vector at a fixed azimuthal angle 
but varies the magnitude with respect to the reciprocal lattice 
vector of the sample. By scanning the scattering vector through 
the Bragg condition in both azimuthal and radial directions we 
can map out the in-plane shape of the Bragg peak. 

In previous SANS experiments on vortex matter, the scatter- 
ing intensities along the azimuthal and radial directions are usu- 
ally resolution-limited. The trade-off for the high resolution of 
a reflectometer is a relatively low neutron flux compared to typ- 
ical SANS instruments. The flux restricts our study to the struc- 
ture of the vortex matter in the low field region of the phase 
diagram where the scattering is strong. 

Additionally, in the typical SANS geometry, one illuminates 
as much of the sample as possible to increase the scattering 
intensity. The measurement is an averaged signal over the bulk 
sample. One particularly relevant aspect of the AND/R geom- 
etry is the width of the beam compared to the sample diame- 
ter (here < 1.0 mm vs. 12.1 mm). Since only a small slice of 



the sample is exposed to the neutron beam, we can shift the 
position of the sample relative to the neutron beam along the x 
direction ( labeled in Fig. 1(b)) and study the vortex matter at 
different sample locations (Hanson et al. 2010). The position of 
the neutron beam relative to the sample (x value) is located by 
scanning through the neutron absorbers (Cd masks) around the 
sample holder and the uncertainty of the x value is ±1.0 mm. 
The data in Fig. 2(a) is a three-dimensional plot of the inte- 
grated intensity versus azimuthal angle cj versus x value, where 
the origin of x represents the center of the sample. The intensi- 
ties have been scaled for the different scattering volumes. From 
Fig. 2(a), it is clear that the Bragg peaks are spatially inhomo- 
geneous (physical implications are discussed in Hanson et al., 
2010). 

4. Three Dimensional Structure Mapping 

The ability to scan the sample along the x direction (see Fig. 
1(b)) allows us to map the Bragg peak from the sample edge to 
the sample center and study the spatial dependence of the vortex 
matter structure. Since the detector is sensitive to the impact of 
neutrons in two dimensions, by adding vertical collimation to 
the instrument, we can observe the orientational disorder along 
the length of the flux lines. 

On top of the standard colUmation, a vertical slit is placed in 
position as shown in Fig. 1(c). The addition of this single verti- 
cal slit in front of the sample focuses our ribbon-shaped neutron 
beam into a point source, located at the slit itself. The distri- 
bution of the neutrons from the source is not completely uni- 
form along the y-direction. This vertical variation is due to the 
fact that the monochromator of AND/R, which is a wavelength 
selector, is composed of seven separate pyrolytic graphite ele- 
ments arranged vertically (Dura et al., 2006). The gaps between 
these elements are visible in Fig. 2(b), a color plot of the inten- 
sity from the main neutron beam on the detector with the ver- 
tical slit in place. Here only the scattering from three fingers is 
measurable with the given aperture widths and positions. 

Fig. 2(c) is the raw scattering intensity at x = —0.75 mm with 
a vertical slit of 1 mm width. In contrast. Fig. 2(d) is the raw 
scattering data from the same state with the vertical sUt wide 
open (20 mm). Comparing these two plots, it is clear that the 
inhomogeneity of the Bragg peak is revealed by the tight verti- 
cal sUt. 

In order to better visualize the orientational disorder along 
the length of the flux lines, we decompose the scattering inten- 
sity on the areal detector into several sections along the y axis. 
The number of resolvable sections depends on the vertical col- 
limation. By back-projecting all possible paths of the neutrons 
from a given j-pixel on the detector to the source, we can map 
the detector regions onto the corresponding vertical locations 
in the sample (z value). In order to correct for the non-uniform 
distribution of the main beam, each section is divided by the 
corresponding main beam intensity (IIImb)- Fig 2(e) shows the 
(corrected) integrated intensity versus vertical position in the 
sample (z = mm corresponds to the bottom of the sample) 
versus the rocking angle cj. In this way we map out both the 
horizontal and vertical variations of the Bragg peak shape. Fig. 
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2(a) and Fig. 2(e) demonstrate that the vortex lattice is spatially 
defective throughout the entire the sample. 

5. Summary and Conclusions 

The high resolution and ribbon-beam geometry of AND/R have 
allowed us to probe the spatial dependence of the vortex matter 
structure transverse to the flux lines in a Nb single crystal. By 
adding vertical pinhole coUimation, we are able to simultane- 
ously measure the orientational disorder along the length of the 
flux hues in addition to the scattering along the azimuthal and 
radial directions. The neutron diffractometer provides us with 
a new approach to map out the complete structure of a vortex 
lattice and enables us to locate the optimum position inside the 
sample in searching for the Bragg glass phase. Spatial variations 
of the orientational order are observed in the lattice structure 
throughout this Nb single crystal. The technique could easily 
be adapted for the spatial mapping of other three dimensional 
structures with suitable geometry. 
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Figure 1 

(a) A sketch of our bulk Nb single crystal in a magnetic field with the corre- 
sponding stmcture of an ideal vortex lattice in reciprocal space, (b-c) The scat- 
tering geometry for the neutron diffractometer. The magnetic field is applied 
parallel to the cylindrical axis (the (111) direction of the Nb crystal), (b) The 
sample and magnetic field are rotated together by an angle ui with respect to 
the incident neutron beam. By using horizontal slits, HS, the incident neutron 
beam is collimated transverse to the flux lines. By shifting the sample along 
the -Y direction, different slices of the sample are probed. An areal detector is 
placed at the angle 28 which is detemiined by the appHed magnetic field. The 
slice in the Nb crystal shows the part of the sample exposed to neutrons, (c) By 
adding vertical collimation, VS, after the horizontal slit in the neutron path in 
conjunction with an areal detector, we can physically distinguish the scattered 
neutrons from different vertical sections of the sample. 




Figure 2 

(a) The integrated intensity of the scattered neutrons from different horizontal 
slicing sections of the vortex matter as a function of ui. The integrated intensity 
is calculated by summing over the areal detector and it has been scaled for the 
different scattering volumes. The x value represents the distance to the center 
of the sample, (b) A two dimensional color plot (on a logarithmic scale) of the 
main beam (26 = 0°) with a vertical slit of 1mm, revealing the dependence 
of the main beam on the monochromator (c) A two dimensional color plot of 
the Bragg peak (28 = 0.23°) aS. x = —0.75 mm with a vertical slit of 1 mm, 
revealing the raw vertical scattering intensity from the vortex matter (d) A two 
dimensional color plot of the Bragg peak (28 = 0.23°) at x = —0.75 mm with 
a wide open vertical slit (20 mm), (e) The integrated intensity from separating 
(c) into seven sections. Due to the vertical variations in the main beam, the inte- 
grated intensity is normalized to the main beam (IIImb)- The z value represents 
the vertical location in the sample. The error bars in (a) and (e) indicate one 
standard deviation from counting statistics. 



